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ABSTRACT 

Most low-mass protostellar disks evolve in clustered environments where they are af- 
fected by external radiation fields, while others evolve in more isolated star-forming 
regions. Assuming that the magneto-rotational instability (MRI) is the main source of 
viscosity, we calculate the size of a poorly ionized, MRI inactive, and hence low viscos- 
ity region - the "dead zone" - in these protostellar disks. We include disk ionization 
by X-rays, cosmic rays, radioactive elements and thermal collisions, recombination by 
molecules, metals, and grains, as well as the effect of turbulence stimulation in the 
dead zone by the active layers lying above it. We also calculate the gap-opening masses 
of planets, which are determined by a disk's viscosity and a disk aspect ratio, for disks 
in these environments and compare them with each other. 

We find that the dead zone is a robust feature of the protostellar disks that 
is largely independent of their environment, typically stretching out to ~ 15 AU. 
We analyze the possible effects of dead zones on planet formation, migration, and 
eccentricity evolution. We show that the gap-opening mass inside the dead zone is 
expected to be of the order of terrestrial and ice giant mass planets while that outside 
the dead zone is Jovian or super-Jovian mass planets, largely independent of the star- 
forming environment. We show that dead zones can significantly slow down both type 
I and type II planetary migration due to their lower viscosity. We also find that the 
growth of eccentricity of massive extrasolar planets is particularly favorable through 
the planet-disk interaction inside the dead zones due to the large gaps expected to be 
opened by planets. 

Key words: accretion, accretion disks - turbulence - planetary systems: formation - 
planetary systems: protoplanetary disks - planets and satellites: general - Solar system: 
formation - stars: pre-main-sequence 



1 INTRODUCTION 

The discovery of over 160 extrasolar planets has opened up 
several fundamental questions about planet formation. Plan- 
etary masses vary from 7.5 Earth masses all the way up to 
~ 15 Jupiter masses. Massive planets are found close to 
their central stars which implies that their migration from 
their points of origin further out in the disk needs to be un- 
derstood. The high eccentricity of many of the planets also 
demands explanation. Theoretical attention has increasingly 
focused on the interaction of planets and their surrounding 
disks in order to elucidate these questions. In this regard, 
while most models still assume that disks possess a constant 
turbulent viscosity throughout, the possibility that "dead 
zones" of very low turbul ence exist in disks has profou nd ef- 
fects on planet formation iMatsumura fc Pudritzl2003l here- 
after MP03). 

* E-mail: soko@physics. monaster, ca (SM); pu- 
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The idea of a dea d zone in protostellar disks was first 
proposed by iGammiel l|l996ft . He argued that there is a 
poorly ionized region where the growth of the magneto- 
rotational instability (MRI) against Ohmic dissipation can- 
not be sustained. Since the MRI is thought to be the most 
promising source of the disk's turbulent "viscosity" , the dead 
zone is expected to have nonzero, but very low viscosity. 

Dead zones have several important effects on planet for- 
mation including (1) planetary masses, (2) planetary migra- 
tion, (3) planet formation via gravitational instability, and 
(4) planetary eccentricity. We discuss these below. 

The first point was ana lyzed in our previous paper 
jMatsumura fe Pudritzl l2005l hereafter MP05): the gap- 
opening mass of planets is expected to be smaller inside the 
dead zone. A planet opens a gap when the angular momen- 
tum transfer rate by the planetary tidal torque exceeds tha t 
of the disk's viscous torque (e.g. iLin fe Papaloizoulll993F) . 
Since a dead zone has a very low viscosity, the gap-clearing 
timescale is sufficiently short compared to the gas-accretion 
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timescale in that region so that a planet won't accrete much 
gas. 

Regarding the second point, dead zones are likely to 
prolo ng the planetary migration timescale (e.g. lChiang et aD 
2002, MP05). This is potentially important because the mi- 
gration timescale estimated in a protostellar disk without 
a dead zone tends to be much shorter than a disk's life- 
time, and as yet we have no obvious mechanism to stop 
planet migration. Protoplanets migrate as a result of their 
resonant interaction with protostellar disks. They tend to 
migrate inward because the magnitude o f the outer t orque 
is usually larger than the inner torque <Wardlll997ft . The 
current theories suggest that planetary migration has essen- 
tially two stages. The first stage is when protoplanets are 
not massive enough to open a gap in the disk; they migrate 
through the disk (type I) as they accumulate mass. The sec- 
ond stage starts when protoplanets become sufficiently mas- 
sive to open a gap; they migrate at the viscous timescale of 
the disk (type II) without accumulating much mass. A dead 
zone can affect both types of migration. A dead zone's effect 
on type II migration can be readily understood - a planet's 
migration speed is significantly slowed down as soon as they 
enter a dead zone. This is because the speed of type II migra- 
tion is directly proportional to a disk's viscosity and because 
the viscosity in a dead zone is expected to be very low. 

A dead zone's effect on type I migration is rather indi- 
rect - planet migration can be halted or even reversed by 
the mass accum ulation at the edge of the dead zone (e.g. 
lThommesf 2005). Even when planets are not stopped at the 
edge of the dead zone, their migration speed is still likely to 
be reduced as soon as they enter it. This is because these 
planets will probably open gaps in the dead zone where the 
gap-opening mass is roughly a couple of orders of magnitude 
smaller. These planets then switch to the type II migration 
that is a f ew orders of magnitude slower than the type I 
migration l)Wardlll99^ . Thus, the dead zone may rescue 
planets from plunging into the central star due to its lower 
viscosity, and may also work as a switch from type I to type 
II migrations or even as a wall to halt or reverse the type I 
migration. We will further discuss this point in §5. 

Regarding the third point, the mass accumulation at the 
edge of a dead zone, or even at the edge of a gap opened by 
a planet, may lead to planet formation via gravi tational in- 
stability, which req uires a dense, cold region (e.g. Bossl ll997l : 
iMaver et alJl20f)3l . The former type of mass accumulation 
occurs due to different evolution speeds between dead zones 
and active zones, while the latter mass accumulation hap- 
pens because the gap-opening timescale tend s to be faster 
than the disk's viscous evolution timescale. iLufkin et alJ 
( 2004) showed that higher density spiral arms formed by 
a massive planet can lead to a subsequent planet formation 
by gravitational instability. 

Regarding the fourth and final point, we propose in this 
paper that the dead zone may be an optimum place for 
the eccentricity evolution, es pecially for a disk - plane t inter- 
action scenario proposed bv IGoldreich fc Hi3 Hq03) . The 
lower viscosity inside the dead zone would make a wider gap 
around a planet, and therefore allow more space and time 
for eccentricity growth. We discuss this point in §6. 

In this paper, we generalize our work in MP05 and cal- 
culate the size of dead zones and planetary gap-opening 
masses in regions of clustered and isolated star-forming en- 



vironments. This allows us to examine how much the size 
of dead zones is changed by its environment, and therefore 
which parameters affect them most. We also discuss each 
possible effect of a dead zone proposed above. We find that 
both the size of a dead zone (typically ~ 10 — 20 AU) and 
gap-opening masses are almost independent of a disk's envi- 
ronment. This is because the disk structures and hence the 
ionization structures in dense regions of disks (e.g. a pos- 
sible region of a dead zone) are almost the same for both 
isolated and clustered star-forming environments. We give 
an overview of a gap-opening mass in §2, and the MRI tur- 
bulence and our disk models in §3. We present our results in 
§4 and discuss dead zone effects on migration and eccentric- 
ity growth in §5 and §6 respectively. Finally, we summarize 
our work in 87. 



2 GAP-OPENING MASS 

The mechanisms of angular momentum transfer in proto- 
stellar disks are not fully understood, but planet formation 
is thought to be linked to at least two of them - tidal interac- 
tion between a gaseous disk and a protoplanet, and viscous 
diffusion of a disk. 

Tidal torques carry away angular momentum most 
efficiently around Lindblad resonances located at about 
a disk pressure scale hei ght h from a protoplanet 
jGoldreich fc Tremainei ri980) and transfer that to the disk 
when generated density waves shock and damp. Assuming 
the disk is inviscid so that there is no viscous torque, and 
that the density waves sh ock immediately in the v icinity 
of the Lindblad resonances, iLin fc Papaloizoul Jl993t) deter- 
mined a gap-opening mass of the planet: 
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(1) 



where M v and M* are mass of a planet and a star respec- 
tively, r v is the orbital radius of th e planet . and h v is the 
pressure scale height at that radius. Rafikovj (120021) argued 
that their assumption is too radical and determined the gap- 
opening mass in an inviscid disk by considering the nonlinear 
evolution of density waves as well as taking account of the 
planet migration effect: 
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where Q is the Toomre parameter. The first term represents 
the case of feedback being not strong enough to stop migra- 
tion - the density waves reflected from the edge of a forming 
gap are not enough to fill in the gap and hence smooth out 
the difference between outer and inner planetary torques. 
The second term corresponds to the case of feedback being 
sufficient to stop migration. All the disk models we used 
turn out to be gravitationally stable (Q > 1), thus planets 
are expected to form through the core accretion rather than 
the gravitational instability, and hence have smaller gap- 
opening masses compared to Equation Q . Equation Q re- 
places Equation @ when density waves damp immediately 
at the Lindblad resonances (at a distance of ~ hp from the 
planet) where they are excited. 

When the disk is viscous, we need to take account of 
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the effect of viscous torque whose primary source is thought 
to be the MRI turbulence. IfBalbusfc Hawlevl <199lft showed 
that a slightly perturbed, weak magnetic field can be ampli- 
fied sufficiently to radially transfe r the angular momen tum. 
Recent numerical simulations fe.g. lFleming et al . 2000) have 
shown that the MRI is active when the dissipation timescale 
of the turbulence (t<jis = h 2 /rf) is more than 10 3 — 10 4 times 
slower than the growth timescale (t grotvt h — Ii/Va)'- 

ReM = M > io 3 _ 10 4 , (3) 
V 

where Va is the Alfven speed and rj is the diffusivity of the 
magnetic field. 1 

In a viscous disk, the gap-opening mass is reached when 
the angular momentum transfer rate b y the tidal torque 
excee ds that by the viscous torque (e.g. iLin fc Papaloizoul 
Il993f> : 

wher e a ss is the viscous a parameter IShakura fc Sunvaevl 
Il973h . For a less viscous disk, or equivalently, a disk with a 
smaller a ss , a planetary gap-opening mass becomes smaller. 
Assuming that the gap-opening mass inside the dead zone is 
well approximated by the gap-opening equation @, MP05 
showed that the gap-opening mass ratio inside the dead zone 
to outside it could be up to about fOO - close to the mass 
ratio of Jovian to terrestrial planets in our solar system. 

Although numerical simulations show th at there is 
some residual gas accretion through a gap (e.g. lKlevlll999l 
iLubow et all H 999). the gap-opening mass i s still expected 
to be close to a planet's final mass. iKlevI lll999l) showed 
that the gas accretion through a gap is markedly reduced 
for a small viscosity: a SB ^ 5 x 10~ 4 . This corresponds to a 
visc osity inside a dead zone found in numerical simulations 
(e.g. iFleming fc Stondl2003h . Thus, we can expect that the 
gap-opening mass is most likely to be the final planetary 
mass inside a dead zone . In a standard disk with mass of 
~ 0.02M^. iLubow et al.l Jl99Sf) showed that the gas accre- 
tion rate through a gap is 10" 5 -10" 6 Mj/year at 10-40 AU 

- it takes 10 5 — 10 6 years to double a mass of Jupiter outside 
a dead zone. This may be shorter than the disk lifetime, but 
longer than the planetary migration timescale in this region 

- type I migration timescale for Jupiter is 10 3 — 10 4 years and 
type II migration timescale for a ss = 0.01 is 10 4 — 10 5 years 
(see Fig. 1101 . Thus, as long as planets migrate inward, they 
probably enter a dead zone and thereby stop gas accretion 
(as they open a clear gap) before they acquire a significant 
amount of mass through gaps. 



1 ISano fc Stonel i2002ft (hereafter SS02) studied the nonlinear 
evolution of the MRI in weakly ionized accretion disks including 
the Hall effect, and obtained the critical magnetic Reynolds num- 
ber of ReM,SS02 = V^/(riQ) ~ 30. Our definition of the magnetic 
Reynolds number relates to theirs as ReM,SS02 = Va/c s ReM = 
1/2 

c*ss Re Mi where c s is the sound speed. Therefore, their critical 
magnetic Reynolds number roughly corresponds to the case of 
ReM = 10 2 for the standard viscosity parameter a aa = 0.01. 



3 MRI AND DISK MODELS 

3.1 MRI — its existence and effectiveness 

Although the MRI is considered to be an effective po- 
tential source of disk viscosity, its existence and effective- 
ness are still under di scussion. Regarding the former point, 
iFromang et alJ B2002T pointed out that a very small frac- 
tion of the cosmic abundance of metal atoms can signifi- 
cantly diminish the size of a dead zone by picking up the 
charges of molecular ions which recombine with electrons 
10 times more effectively than metal ions (i.e. keeping the 
overall ionization rate high). The electron fraction, however , 
can be significantly reduced by grains (e.g. lSano et ai]|2000l) 
because they recombine with electrons even more efficiently 
than molecular ions. This effect can potentially make the 
midplane region of an entire protostellar disk MRI dead. We 
consider recombination rates of electrons with metal ions, 
molecular ions as well as grains as in MP05 (also see Ap- 
pendix A in this paper). To calculate the ionization rates, 
we take account of X-rays, cosmic rays, radioactive elements 
and thermal collisions of alkali ions (for details, see MP03 
& MP05). 

Dead zones have been studied by many authors 

both analytical ly (e. g. iGammid 1996; pano et alJ l2000l : 

iGlassgold et al.l 120001: IFromang et al.l l200l MP 03, MP05) 
and numerically (e.g. IFleming fc Stond l2003tl . Recently, 
llnutsuka fc Sana (|2p05) suggested that once MRI turbu- 
lence becomes active, it can be sustained in a protostel- 
lar disk without any external ionization and concluded that 
most regions in protostellar disks remain magnetically ac- 
tive. For dense regions of a disk with dust grains, they pro- 
posed that collisions of energetic electrons might provide 
enough ionization to sustain MRI. In considering their argu- 
ments however, we note that their estimate for the electron 
fraction is still one to two orders of magnitude smaller than 
the one required by recent numerical work. For a disk sur- 
face, assuming dust grains are depleted from the region, they 
proposed that turbulent eddies could mix the ionized region 
with the neutral region, and therefore homogenize the ion- 
ization region, if the recombination rate is sufficiently low. 
Grains could indeed be absent in a high temperature re- 
gion close to the central star (see §4.3). However, both their 
effective electron recombination coefficient and the initial 
electron fraction are one to two orders of magnitude too 
small because of their assumption that a as = 0.2. 2 Thus, 
we consider these processes are probably unable to sustain 
the MRI turbulence throughout the entire disk. 

Even if the MRI is present, other effects like ambipolar 



2 Using the formula given in Inutsu ka fc Sancl |2005|> . the effec- 
tive recombination coefficient can be written as & ~ 10~ 7 " t+ + 

n e 

10~ 12 , where n e , n m +, and n M+ are number densities of 
electrons, molecular ions, and metal ions respectively (see Ap- 
pendix A). At regions close to the central star, densities of elec- 
trons and metal ions are about two orders of m agnitude larger 
than the density of molecular ions iSano et alj|200ol) . This will 
gi ve 0' ~ 10~ 9 and the el ectron fraction of ~ 10~ 14 (see eq. (12) 
in llnutsuka fc Sandl20oih . unless the initial electron fraction is 
larger. Also, their suggested electron fraction to sustain MRI is 
likely an underestimate because they adopted a relatively large 
viscosity parameter ce BS = 0.2. 
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diffusion could work against its angular momentum trans- 
port. Usin g three-dimensional m agneto-hydrodynamic sim- 
ulations, lHawlev fc Stone! l)l998f) showed that there is sig- 
nificant angular momentum transfer when the ion-neutral 
collision rate is 100 times larger than the local epicyclic fre- 
quency, and that ions and neutrals are essentially decoupled 
when this ratio goes below 0.01. We calculated the corre- 
sponding ratio in our disk models and found that the value 
is typically ~ 10. Therefore, we expect at least some angular 
mo mentum transport i n acti ve regions of our disk models. 

iFleming fc Stone! (|2003ll studied the evolution of the 
MRI in stratified accretion disks and found that there is a 
minimum level of angular momentum transport in the dead 
zone in the presence of active layers due to the Reynolds 
stress. In the case relevant to ours (Ejjz/Eals ~ 10), they 
obtained a ss ~ 10~ 3 for active layers and a ss ~ 5 x 10 -5 
for a dead zone. Later we will show that this a 3S value 
doesn't affect the gap-opening mass inside the dead zone. 
Another effect discussed in the literature is the mass mix- 
ing betweenji_dead zoneand the active layers lying above 
it jFleming fc Stond l2003h . Although this vertical mixing 
may have no effect in radial angular momentum transfer, 
we take account of this effect by imposing the condition 
E_dz/Sals 10 on a dead zone besides ReM < 10 3 — 10 4 . 



3.2 Disk Models - Clustered vs Isolated Star 
Formation 

Typical protostellar systems are expected to form in a clus- 
ter like Orion Nebula, where they are irradiated by nearby 
luminous OB stars. On the other hand, some star formation 
regions li ke Taurus are known not to have young massive 
stars (e.g. iLuhman et alJl2003^> . It is interesting to compare 
these two star formation regions from the viewpoint of their 
planet formation environments. 

To determine the size of dead zones of disks in these 
environments, we need to calculate ReM (see equation (J^J) 
throughout disks that depends both on the thermal struc- 
ture of the disks via h and on the ionization structure via r). 
Both of these quantities will be affected by the presence of 
massive stars. 

The difference in the thermal structure for protostellar 
disks in the clustered and isolated case is that the former 
disks are exposed not only to their central stars, but also to 
other massive stars in the cluster. For the isolated disks, as in 
MP03, we adopted the isolated , two-layered disk models de- 
veloped bv lChiang et alJ J200lf) . For the disks in a cluster, as 
in MP05, we used the disk models bv lRobberto et alJ (2002) 
where they modified the iso lated two-layered disk models of 
IChiang fc Goldreichl il997f) by submerging these disks in a 
cluster environment. In their model, the external star is as- 
sumed to be located at an average distance of 0.1 pc from 
a low-mass protostellar disk and have a stellar luminosity 
of L, — 6 x 10 38 ergs s _1 . Due to this extra heating, the 
disks in a cluster are more flared compared to the isolated 
cases (see Fig. 0. In both disk models, disk temperatures 
are calculated by assuming a radiative, hydrostatic equilib- 
rium. For a particular disk radius, these models give two 
temperatures - one in the surface layer, and the other in the 
disk interior. Also, both disk models assume a power-law 
surface mass density: E = Eo (xjf) 32 ' ■ 



Fig. compares temperatures of these models (left 
panel) as well as corresponding disk aspect ratios (right 
panel). From the left panel, it is apparent that the disk inte- 
rior temperatures (i.e. the mid-plane temperatures) are the 
same for both environments out to ~ 10 AU, while the radi- 
ation from external fields clearly affect on a disk's thermal 
structure beyond ~ 10 AU, giving a nearly constant temper- 
ature (~ 50 K) there. The temperature difference of up to a 
factor of ~ 2 leads to a more flared disk, and hence a larger 
disk aspect ratio and larger gap-opening masses in a clus- 
tered star- forming environment (see Equation (2) and @). 
From the right panel, we can see that, at 20 AU, the disk as- 
pect ratio for isolated and clustered star-forming regions are 
0.09 and 0.15 respectively. This ~ 70 % difference in a disk 
aspect ratio leads to a gap-opening mass of 0.8 and \.2Mj 
(where Mj is a mass of Jupiter) in each environment for a 
standard alpha parameter of a sa — 0.01 (see also: Fig. 0. 
Note that, by knowing the radial temperature structure of 
the disk, and hence its aspect ratio; h/r, the only parameter 
needed to determine the distribution of gap-opening masses 
is a disk's viscosity parameter a ss . We will discuss this in 
§4.1 and 4.2. Also, the temperature at the innermost disk 
reaches ~ 1500 K, suggesting that the region may be void 
of dust grains due to thermal evaporation. We will consider 
this effect in §4.3. 

Note that, although the estimated surface layer's tem- 
perature of the disk model in a clustered environment is 
very low; ~ 200 K beyond 10 AU, compared to the at- 
mosphere temperature assumed for photoevaporating disk 
mode ls; ~ 10 4 K (e.g. IShu et alJ 1 19931: iHollenbach et alJ 
Il994fl . it does not affect our results so much. This is because 
the disk atmosphere reaches that high temperature only in 
very high altitude regions - about one order of magnitude 
higher than our surface disk height, where the density is very 
low. In addition, the gap-opening masses depend on the disk 
interior temperature, which is related to the pressure scale 
height rather than the disk surface temperature. 

The difference in the ionization structure for protostel- 
lar disks in the clustered and isolated case is that the former 
is also exposed to the X-ray radiation from an external star 
besides other ionization sources common in both cases (see 
below). Here, we assume the X-ray luminosity of an exter- 
nal star to b e L x = 10 34 ergs s - 1 , which is likely to be an 
upper limit (IStelzer et al.ll2005T) . 3 In both environments, 
the ionization of disks is due to X-rays fr om the central 
star (L x = 10 30 ergs s _1 with kT x = 2 keV. iFeigelson et alJ 
(2002)), cosmic rays {$,cr ~ 10~ 17 s" 1 with the attenua- 
tion length of ~ 96 g cm~ 2 . ISano et al l i2000|) ), radioactive 
eleme nts (£ha ~ 6.9 x 10~ 23 s" 1 . lUmebwashTfc Nakand 
( 1981)), as well as heated alkali ions (important in the high 
temperature, inner part of the disks). 



3 We performed the simulation with a more typical value of the 
X-ray luminosity of a massive star: L x = 10 31 ergs s _1 , and found 
that the external star has less effect on the total ionization rate. 
But this does not change the extent of the dead zones. 
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4 DEAD ZONES AND GAP-OPENING 
MASSES 

In this section, we compare protostellar disks in different 
environments and obtain sizes of their dead zones (§4.1) as 
well as the gap-opening masses of planets (§4.2). Our fo- 
cuses are on comparing (1) disks in a stellar cluster and 
isolated disks, and (2) disks with and without cosmic ray 
ionization. The latter case is included because low energy 
cosmic rays, which are more responsible f or ionization com- 
pared to higher energy cosmic rays (e.g. lLepdll992T) . may 
be partially excluded fro m the disks by magnet ic scattering 
jSkilling fc Strondll97r3) . Note, however, that iDesch et alJ 
i2004) recently showed that galactic cosmic ray radiation is 
likely to be abundant in protostellar disks. Here we include 
this comparison for completeness. In the end of this section, 
we also discuss the effect of dust evaporation (§4.3). 

4.1 Clustered vs Isolated Star Formation — Dead 
Zones 

Fig. [5] shows dead zones in all representative cases of our 
study - disks in a cluster, isolated disks, and disks with 
dead zones estimated by X-ray ionizations alone for each en- 
vironment. Although we cover a large parameter space: the 
surface mass density of E = 10 s - 10 4 g cm , the magnetic 
Reynolds number of Re M = 10 2 - 10 4 and a aa = 1.0-0.001, 
we only show the case of (ReM, a ss ) = (10 3 , 0.1 — 0.001) 
with Eo = 10 3 g cm -2 in this figure. We will discuss the 
effects of different parameters below (see Fig. [3] as well). 

Comparing disks in a cluster with isolated disks (upper 
two panels), we find that dead zones stretch out to ~ 10 — 20 
AU in both environments. This is not so surprising because 
the effect of an external star and a nebular environment is 
not dominant within ~ 10 AU (see Fig. 0. The insensitiv- 
ity of the dead zone size to thermal environments suggests 
that protostellar disks in both environments have a similar 
transition radius of disk's viscosities. 

Comparing the dead zones estimated by the total ion- 
ization (upper panels in Fig. |5J with those estimated by 
the X-ray ionization (lower panels), we can see that cosmic 
rays have a large effect on a dead zone shape at outer radii. 
Disks ionized by both X-rays and cosmic rays (as well as 
other sources) have fully magnetically active regions beyond 
~ 10 AU, while disks ionized by X-rays alone are expected 
to be able to sustain MRI only when there is a significant 
mass-mixing between active and dead zones (note that the 
dashed line represents a critical surface mass density ratio of 
the dead zone and the active layers below which mass-mixing 
between these two zones is not negligible) . This is essentially 
due to the geometrical difference between these two major 
ionizing sources, X-rays from a central star and cosmic rays: 
X-rays are emitted from a stellar magnetosphere and there- 
fore penetrate disks with an angle of <C 90°. As a result, 
disks become more optically thick toward the outer part of 
a disk for X-rays. Cosmic rays on the other hand, will prop- 
agate down to the disk preferentially along disk's magnetic 
field lines that are orthogonal to the disk's surface, hit the 
disk with an angle of ~ 90° and experience a similar optical 
thickness throughout the disk. 

As described in MP05, we define edges of the dead 
zones as intersections of a critical surface mass density ra- 



tio (dashed line) with curved boundaries of dead zones. The 
dead zone edges obtained this way are plotted in Fig. [3] 
where we can compare parameters' effects on the dead zone 
sizes in different environments and ionization sources. In a 
standard parameter range (ReM, a ss ) = (10 3 — 10 4 , 0.1 — 
0.001), dead zones estimated by the total ionization (upper 
two panels) tend to be about < 10 — 20 % larger for disks 
in isolated environments. For ReM and a 3S , the size of dead 
zones changes by about a factor of 2 each over three orders of 
magnitude. The most significant effect however comes from 
a surface mass density of the disk at 1 AU Eo, for which 
the dead zone's size changes by the same factor, 2, over one 
order of magnitude. 

The difference in a dead zone size due to environ- 
ments becomes more apparent for disks ionized only by X- 
rays (lower two panels). This is because of the extra X- 
ray ionization effect by a nearby massive star in a clus- 
tered star-formation region. In a relatively heavy disk with 
Eo = 10 4 g cm -2 , dead zones of isolated disks are up to 
about a factor of 2 larger than those of disks in a cluster. 
For ReM and a ss , the size of dead zones changes by about 
a factor of 2 and 1.5 respectively over three orders of mag- 
nitude, while for Eo, the dead zone size changes by a factor 
of > 3 over just one order of magnitude. Thus, the envi- 
ronmental effect becomes more important if cosmic rays are 
excluded from the disks. Note however, that our estimated 
X-ray ionization rate for a clustered environment is proba- 
bly overestimated, since we have not taken account of any 
extinction for X-rays traveling to the protostellar disk and 
since the adopted X-ray luminosity for the massive external 
star is likely to be an upper limit as mentioned earlier. 

In all cases, the most important parameter for the dead 
zone size is the surface mass density of the disk at 1 AU Eo. 
Both ReM and a ss also have a large effect on the dead zone 
size, while the environmental effect on the dead zone seems 
to be negligible unless the X-ray ionization from external 
sources is very powerful. 

4.2 Clustered vs Isolated Star Formation — 
Planetary Gap-opening Masses 

The general effect of a dead zone upon planetary masses is 
to create a bimodal distribution. Gap-opening masses will 
be sharply reduced as we move inwards through the well- 
coupled outer region of a disk and encounter the dead zone. 
This is shown in Fig. 2] which plots the gap-opening mass 
as a function of disk radius by assuming a ss = 0.01 in the 
active region. The figure shows the gap-opening mass for the 
case of a ss = 0.01 throughout the disk (upper dashed curve), 
as well as that for the case of inviscid disks (lower dashed 
curve). The actual distribution is the heavy black curve, 
which follows the low, inviscid curve in the dead zone, and 
then precipitously jumps up to the upper curve in the active, 
well-coupled outer region of the disk. The maximum mass 
of a planet in the disk's dead zone is M maXj DZ ~ 0.09Mj in 
this example. Just outside of the dead zone, the minimum 
mass in the well-coupled active zone is M m i n ,Az ~ 0.79Mj. 

We generalize this approach in Fig. [£|which shows gap- 
opening masses as a function of a disk radius for disks in 
a cluster (left panels) and isolated disks (right panels). In 
each panel, we plot the gap-opening mass lines for disks with 
different, constant values of a as (1.0, 0.1, 0.01 and 0.001) as 
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well as Rafikov's inviscid case (see Equation ©). The dead 
zone radius for a disk with each value of a sa depends on 
the Reynolds number ReM (MP05). Thus, a cross appears 
on each a 3S = constant curve, which locates the dead zone 
radius for that particular disk model. Since we show the data 
for these 3 different values of ReM, each panel of this figure 
has 3 nearly vertical lines, which represent the locus of dead 
zone radii for all disk models for each value of ReM. These 
loci show that the dead zone radii are not highly sensitive 
to the precise value of the magnetic Reynolds number ReM- 

The reader may construct the predicted gap-opening 
masses for any disk model (characterized by a value for a ss 
in the well-coupled active zone, magnetic Reynolds number, 
surface mass density, and ionization) by choosing the panels 
in Fig. to construct a figure akin to that of Fig. ^ One 
of the examples of gap-opening masses throughout the disk 
is plotted in a heavy black line in each panel. Also plotted 
is a mass inside a Hill radius (see Equation (0) which is 
probably an upper limit for a gap-opening mass in an invis- 
cid disk. This line also gives a rough estimate of planetary 
masses formed via disk instability. If the disk's viscosity is 
larger than ~ 0.01, these planets are expected to keep on 
accreting gas after their formation. 

An important implication of the results is that Jupiter 
or more massive planets cannot be formed inside a dead zone 
(within ~ 10 — 20 AU) through core accretion, because a 
gap opens for much lower mass planets - even for terrestrial 
planets. These panels also indicate that an inviscid region in 
a typical protostellar disk gives the same gap-opening masses 
as those expected in a viscous disk with a ss = 10 -4 — 10 -5 . 
Therefore, disks with a SB = 10~ 4 — 10~ 5 can be treated as 
dead. 

Another implication of Fig. |S]is that super massive Jo- 
vian planets are difficult to form. For example, to obtain 
10Mj within a reasonable radius ~ 40 AU, the required 
viscosity parameter is a ss =0.1 for a clustered environment 
and ~ 1 for an isolated environment. Therefore, we speculate 
that those massive planets may be formed directly through 
the gravitational instability of the disk, or they could accrete 
much more gas as they migrate. Regarding the latter case, 
a moving planet is shown not to deplete its feeding zone 
as m uch as a static planet (e.g. Rafikov 200^; lAlibert et alJ 
120041) . 

We saw in Fig. 2] that the bimodal character of plane- 
tary masses created by dead zones can be characterized by a 
large jump in planetary masses; M max ,DZ — > M m i n ,Az- Fig.HJ 
shows this jump in gap-opening masses just inside and out- 
side the dead zone radii. For each value of ReM and a ss , the 
dead zone radii seen in Fig. |3] can be combined with gap- 
opening masses in Fig. |S] to produce Fig. |S] These panels 
clearly show there is a jump in mass at the edge of a dead 
zone. Since dead zone sizes are more or less the same in both 
clustered and isolated environments (see Fig. [3J|, and since 
gap-opening masses around typical dead zone radius (a few 
tens of AU) are about the same (see Fig.QjJ, these minimum 
and maximum gap-opening masses at a dead zone radius are 
about the same for both clustered and isolated star-forming 
environments. 

For a heavier disk with Eo = 10 4 g cm -2 (right panels), 
the "minimum" and "maximum" gap-opening mass lines in- 
tersect for a small viscosity parameter of q ss ~ 10 -3 . This 
is because gap-opening masses for an inviscid disk exceeds 



those for a viscous disk for a sufficiently small value of a S3 
(see the intersection of the lower two curves in the lower 
panels of Fig. 15^. Therefore, in a moderately heavy disk with 
a ss ^ 10~ 3 , there is no jump in gap-opening masses. 

Environmental effects on gap-opening masses can also 
be seen in these heavier mass disk cases. In the upper right 
panel of Fig. |5] where the dead zones are calculated by the 
total ionization, gap-opening masses are significantly differ- 
ent compared to the corresponding case for Eo = 10 3 g cm -2 
(upper left panel), despite the fact that the dead zone sizes 
are about the same in both environments in both cases. This 
implies the difference in gap-opening masses between two en- 
vironments seem more significant in a heavier disk. This is 
because dead zone radii for Eo = 10 4 g cm -2 is about 20—50 
AU, while those for E = 10 3 g cm" 2 is about 10 - 25 AU, 
and because environmental effects kick in only beyond ~ 10 
AU. Another example is a heavy disk ionized only by X-rays 
(lower right panel). Here, gap-opening masses in two envi- 
ronments give more or less the same values. However, note 
that we are comparing these masses at very different dead 
zone radii - for example, for (ReM, ct ss ) = (10 3 , 0.01), the 
dead zone radius is 36 AU for a clustered environment and 
61 AU for an isolated environment. Therefore, these results 
imply that gap-opening masses in an isolated environment 
is changing less sharply, as we can see in Fig. |S] 

In the left panel of Fig. [7| we compare gap-opening 
masses at 20 AU in clustered and isolated star-forming 
regions. The gap-opening masses in a clustered environ- 
ment are about a factor of 1.5 larger than the ones in 
an isolated environment. For standard viscosity parame- 
ters ct S s = 10 -3 — 0.1, the gap-opening masses change from 
roughly a Saturn mass to a few Jupiter masses at 20 AU. 
This indicates the importance of the viscosity parameter on 
a gap-opening mass. 

In the right panel of Fig. |7| we plot the gap-opening 
mass ratios of clustered and isolated environments for both 
viscous (heavy line, see Equation Q) and inviscid (light 
line, see Equation J5J) cases. Note that the result for a vis- 
cous case is independent of an actual value of a viscosity 
parameter. Different star-forming regions change a plane- 
tary gap-opening mass by a factor of ~ 1.2 — 2.8 beyond 10 
AU in both viscous and inviscid disks. This result further 
confirms that the environmental difference is probably not 
a large factor in determining gap-opening masses. 



4.3 Grain Evaporation in the Inner Disk 

Until now, we have neglected the effect of grain evaporation 
at a high temperature region within the innermost region 
of an accretion disk. We take account of grains' thermal 
evaporation temperature (T = 1750 K for graphite g rains, 
T = 1400 K for silicate grains: Iflillenbrand et alJll992l) . and 
calculate the size of a dead zone for ReM = 10 3 cases. The 
dead zones now look as in Fig.|S| We have an inner active re- 
gion due to thermal collisions of alkali ions besides the lack of 
grains in that region. The corresponding gap-opening masses 
are plotted in Fig. [5] Note that in situ planet formation in 
the inner active region (r ^ 0.04 AU) is unlikely. This is 
because the region has no dust (therefore core accretion sce- 
nario is impossible) and is gravitationally stable (therefore 
disk instability scenario is impossible). This inner turbulent 
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region is likely to accrete onto the star rather quickly and 
form an inner hole due to a higher viscosity. 



5 DEAD ZONES AND PLANETARY 
MIGRATION 

In §1, we argued that the dead zones can slow down plane- 
tary migration. Since the type I migration does not depend 
on a disk's viscosity, the dead zones cannot directly affect 
the low mass planet migration. However, when a low mass 
planet starts its migration from outside the dead zone, it can 
be strongly affected by the existence of the dead zone. This 
is because of the difference in disk evolution speed between 
a dead and an active zone. Due to the smaller viscosity in a 
dead zone, the disk mass accreting toward the central star 
from an active zo ne is likely to b e accumulated at the edge 
of the dead zone iGammielll996l) . If a migrating planet sees 
this denser region, the inner torque could be comparable to, 
or even grater than, the outer torque in magnitude. This 
will stall the planet, or even reverse its migration. As we 
argued in §1, the type I migration of a non-gap-opener can 
in principle be divided into two groups: (1) when there is 
enough mass accumulation at the edge of a dead zone, a 
protoplanetary migration is either stalled or reversed, and 
(2) when there is not enough mass accumulation, a planet 
will migrate into a dead zone and likely to open a gap. 

Roughly speaking, for a planet to enter a dead zone 
without being stopped by the mass accumulation at the 
dead zone edge, it has to migrate fast enough compared to 
the disk's mass at around the edge of a dead zone. Fig. I1UI 
compares type I migration timescales of Jupiter, Saturn, 10 
Earth mass planet, and Earth (heavy lines, from bottom to 
top) with disk's viscous evolution timescales (or equivalently 
type II migration timescales) for a SB = 1 — 1CP 5 (light lines, 
from bott om to top). Her e we used the following timescale 
equations jTeraueml l2003) 




where Me is the Earth mass. For a S3 = 0.01, it is apparent 
that Saturn or more massive planets migrate much faster 
than disk's mass while 10 Earth mass or lighter planets 
migrate slower. We expect that this mechanism may keep 
lighter planets beyond a dead zone and help them to grow 
larger. 

A dead zone's effect on type II migration can be read- 
ily seen in Fig. 1101 For example at 10 AU, the migration 
timescale is > 10 4 years for a BS = 0.01 (in an active zone), 
and > 10 6 years for a ss = 10~ 4 (in a dead zone). Thus, type 
II planets slow down significantly as soon as they enter the 
dead zone due to the lower viscosity there. 



6 DEAD ZONES AND PLANETARY 
ECCENTRICITY 

One of the major surprises about extrasolar planetary sys- 
tems is the presence of isolated planets with large orbital 



eccentricities. Moreover, all very massive extrasolar planets 
(M p > 7Mj) found so far have rather large eccentricities 
(e > 0.2). We argue that this phenomenon may also be ex- 
plained by the existence of a dead zone. 

Current models suggest that the growth of plane- 
tary eccentricity takes place while the disk is still present. 
In such a situation, eccentricity evolution can occur ei- 
ther through planet-planet interac tion (e.g. |CJiiarrg et_alJ 
|2002j) or planet-d i sk int eraction (e.g. lGoldreich fc Sari! 2003: 
bari fc Goldreichl 12004 hereafter GS03 and SG04 respec- 
tively). Planet-planet interaction would require similar mass 
(or more massive) planets to enhance their eccentricities, 
while all systems with a very massive planet do not have 
such a companion. It is possible that a similar mass com- 
panion was ejected out of the system due to a strong in- 
teraction. The other possibility - planet-disk interaction - 
can increase the eccentricity if a planet with an initial mild 
eccentricity opens a large enough gap (SG04). Since larger 
gaps are expected in regions of smaller viscosity, we suggest 
that planets caught in gaps within dead zones are likely to 
grow their eccentricity. 

Conditions during gap formation may be particularly 
favorable for eccentricity growth of a planet (GS03, SG04). 
Planetary interaction with disks at Lindblad resonances en- 
hances their eccentricity whereas interaction at corotation 
resonances damp it. A necessary condition for eccentricity 
growth therefore is that the non-linear saturation of the 
corotation resonances occur. This condition implies that the 
timescale on which the density gradient is flattened by coro- 
tation resonance (t sa t), is shorter than both the disk viscous 
timesclale (t V is) and the timescale on which the gap of width 
w is opened by the principal Lindblad resonances (t gap ); 
tsat < (tvis, tgap) (SG04). It is also necessary that the eccen- 
tricity must grow faster than the gap grows. Slightly mod- 
ifying their paper, we show that the initial eccentricity has 
to be at least 

eo.min = max [ei jm i n , e 2 , min] (7) 



e2,min = 




for this kind of eccentricity evolution to occur (see Appendix 
B for the derivation). 

Fig. 1111 shows a calculated minimum initial eccentric- 
ity for Jupiter and 10 Jupiter mass planet in the disk of 
Fig. [SJwith a ss — 0.01. It is immediately clear that the re- 
quired initial eccentricity for very massive planets is very 
low (eo.min ~ 10~ 3 ) inside a dead zone while that outside a 
dead zone is more than 0.1. This suggests that eccentricity 
of heavy planets may have been excited inside dead zones 
through disk-planet interaction as they migrated into them 
and opened gaps. It is interesting to note that all high ec- 
centricity (e > 0.2) planets observed so far, including very 
massive planets (Mp > 7Mj), are found between 0.07 — 10 
AU - roughly the region of a dead zone. If a planet's ec- 
centricity is enhanced due to an interaction with a disk, we 
expect planets with small eccentricities to be found beyond 
10 AU. 

The condition that ei, m j n = e2, m in in Equation J7J gives 
the absolute minimum of the initial eccentricity required to 
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excite eccentricity growth. This occurs for roughly a Uranus- 
mass planet. We can determine the maximum and minimum 
mass planets whose initially required eccentricity is reason- 
ably small (e.g. eo = 0.01). By setting ei jIn i n = 0.01, and 
assuming a disk mass of Md = 0.01M© as well as a dead 
zone's viscosity of a ss — 5 x 10 -5 , we can determine the 
maximum mass planet whose eccentricity is likely to be ex- 
cited through planet-disk interaction. We find that all very 
massive planets (up to ~ 20Mj) may grow their eccentric- 
ity inside the dead zone in this way. Similarly, by setting 
e2,min = 0.01, we can determine the minimum mass planet. 
Our calculation shows that planets less than an Earth-mass 
in a disk of Md = O.OIM© are too light to interact with a disk 
and enhance their orbital eccentricity. Combining these re- 
sults together, we can see that planets with masses between 
Earth and 20 Jupiter masses may enhance their eccentricity 
reasonably easily inside the dead zone. 



disk interaction. As we discussed in §6, this may be espe- 
cially true for very massive planets (M p > 7Mj). 

In our future work, we revisit the problem of planetary mi- 
gration and dead zones by using time-dependent calcula- 
tions. 
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7 CONCLUSIONS 

We calculated the sizes of dead zones and gap-opening 
masses in models of stationary protostellar disks. We 
performed an extensive parameter search: Eo = 10 3 — 
10 4 g cm -2 , a sa = 10~ 3 — 1, Re M = 10 2 — 10 4 , and compared 

(1) disks in a clustered environment and isolated disks, and 

(2) disks with and without cosmic ray ionization. Our major 
findings are as follows: 

(i) Dead zones are robust features of protostellar disks. 
We took account of all the major sources of ionization and 
recombination and considered other effects like ambipolar 
diffusion and mass mixing between active layers and a dead 
zone. We showed that dead zones typically stretch out to a 
few tens of AU. 

(ii) Cosmic ray ionization has a large effect on an ioniza- 
tion structure at outer part of the disk (beyond 10 AU), but 
does not change the sizes of dead zones significantly. From 
Fig. |21 it is clear that the ionization of an outer part of the 
disk is dominated by cosmic rays. The sizes of dead zones 
determined by the total ionization and those by the X-ray 
ionization are about the same for the same parameter set 
(So, a sa , ReAf), in the same environment (clustered or iso- 
lated star- forming regions) , because the ionization structure 
in inner dense regions are similar (i.e. dead zone boundaries 
intersect with a critical mass ratio at similar radii) . The dif- 
ference due to ionization sources become more apparent in a 
denser disk (Eo ^ 10 4 g cm -2 ) because disks become more 
optically thick. 

(iii) The size of a dead zone depends mildly on magnetic 
features of a disk (a sa and Rem), and rather strongly on a 
disk surface mass density Eo (see Fig. |SJ . It is almost inde- 
pendent of disk's environments (clustered or isolated star- 
forming regions). 

(iv) Jovian or super Jovian planets are likely to be formed 
beyond a dead zone. Inside dead zones, a gap opens for 
smaller mass planets - ice giants or even terrestrial plan- 
ets (see Fig. 01 and . 

(v) Dead zones may significantly slow down both type I 
and type II planet migration (see §5). 

(vi) Gaps within the dead zones may be good regions in 
which planets enhance their orbital eccentricities via planet- 
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Figure 1. Left: temperature structures of disk models in clustered (heavy lines) and isolated (light lines) environments. Solid lines show 
disk temperatures below the pressure scale height of the disk, and dashed lines show disk envelope's temperatures. Right: disk aspect ratios 
of disk models in clustered (heavy line) and isolated (light line) environments. External fields' effects become important in an outer part 
of the disks. 
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Figure 2. Comparing dead zones in clustered (left panels) and isolated (right panels) environments. For the upper panels, ionization sources 
are X-rays from the magnetosphere of a central star and an external star (for a disk in a clustered environment) as well as cosmic rays, 
radioactive elements, and thermal ionization by alkali ions. For the lower panels, the ionization sources are X-rays from the magnetosphere 
of a central star and an external star (for a disk in a clustered environment). Here, we choose Eo = 10 g cm" " 2 with Re M = 10 3 . The 
uppermost line shows the surface disk height; the lowermost line shows the pressure scale height, while three curves indicate the dead zone 
boundaries for a ss = 0.1, 0.01, and 0.001 from inside to outside. The dashed line is where Sbelow/^above = 10- 



12 Soko Matsumura and Ralph E. Pudritz 

e = 10 3 g cm- 2 





isolated 


-- Re M =10 4 


clustered 


■■■ He M =10 3 


isolated 


■■■■ Re M =10 3 


dllHtlTCd 


Rc M =10- 


isolated 


— Re M =1 ° 2 


clustered 



Viscosity parameter a 



5 

< 



Re M =10 , isolated 
Re H =10 .clustered 

Re M =10", isolated 
■ ■ ■ ■ Re M =10" , clustered 
— Re M =10\ isolated 
Re = 1 0~, clustered 




Viscosity parameter a 



Re M =IO . isolated 
Re M =IO , clustered 
Re M =lt) 3 , isolated 
Re M =10", clustered 
Re M =10~, isolated 
Re =10", clustered 



Viscosity parameter a. 



Total 



Re M =IO . isolated 
Re M =10 , clustered 
Re M =10 3 , isolated 
Re M =10 3 , clustered 
Re M =10 2 , isolated 
Re.=10~, clustered 




X-rays 



Viscosity parameter a 



Figure 3. Comparing dead zone radii in clustered (heavy lines) and isolated (light lines) environments. Solid, dotted and dashed lines are 
for ReM — 10 2 , 10 3 , and 10 4 respectively. The difference between two environments are more apparent in the cases without cosmic ray 
ionization (lower panels) compared to general cases (upper panels). Note the difference in scale of the y axes. 
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Figure 4. Gap-opening masses for a disk exposed to an ex- 
ternal star which has the surface mass density, at 1 AU, of 
So = 10 3 g cm~ 2 . The lower dashed line shows the gap-opening 
mass for the case of no MRI viscosity, while the upper dashed line 
shows the gap-opening masses for a ss = 0.01. For the magnetic 
Reynolds number Rcm = 10 3 , the fiducial minimum gap-opening 
mass throughout the disk is shown in a heavy line. 
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Figure 5. Comparing dead zones in clustered (left panels) and isolated (right panels) environments. Dead zone radii are calculated by 
using the total ionization rate in all cases. For upper panels, we assume the surface mass density at 1 AU of Sq = 10 3 g cm -2 , while for 
lower panels, So = 10 4 g cm -2 . The lowermost line shows the gap-opening mass for the case of no MRI viscosity, while the upper parallel 
lines show the gap-opening masses for a different strength of magnetic field: the solid line is for a sa = 1, the dotted line is for 0.1, the 
dashed line is for 0.01, and the dot-dashed line is for 0.001. The dead zone for the magnetic Reynolds number, ReM = 10 2 , is inside a solid 
line with crosses, that for 10 3 is inside a dotted line with crosses, and that for 10 4 is inside a dashed line with crosses. The thick solid line 
traces our fiducial minimum gap-opening mass throughout the disk for MRI "viscosity," a ss = 0.01. 
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Figure 6. Comparing gap-opening masses just inside and outside the dead zone radii in clustered (heavy lines) and isolated (light lines) 
environments. Solid, dotted and dashed lines are for Re^ = 10 2 , 10 3 , and 10 4 respectively. The difference between two environments are 
more apparent in the cases without cosmic ray ionization (lower panels) compared to general cases (upper panels). Note the difference in 
scale of the y axes. 
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Figure 7. Left: comparison of gap-opening masses at 20 AU in clustered (heavy line) and isolated (light line) environments. Note that the 
magnetic Reynolds number does not affect the gap-opening masses. Right: gap-opening mass ratios of isolated and clustered star-forming 
environments. Heavy line is for a viscous disk (see Equation JIJ) and light line is for an inviscid disk (see Equation 0). Note that the 
results for a viscous disk are independent of a viscosity parameter ct SB . 
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Figure 8. The effect of grain evaporation on the sizes of dead 
zones calculated for a protostellar disk exposed to an external 
star. The ionization sources are X-rays from the magnetosphere 
of a central star and an external star as well as cosmic rays, 
radioactive elements, and thermal ionization by alkali ions. We 
use So = 10 3 g cm -2 and Rem = 10 3 with a ss = 0.1, 0.01, and 
0.001. For the explanation of each line, see Fig. l2l 
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Figure 9. The effect of grain evaporation on gap-opening masses 
for a disk exposed to an external star which has the surface mass 
density, at 1 AU, of S = 10 3 g cm . For the explanation of 
each line, see Fig. l5l 



Figure 10. Migration timescales as a function of disk radius. 
Heavy lines show, from bottom to top, type I migration timescales 
for Jupiter, Saturn, 10 Earth mass planet, and Earth respectively. 
Light lines show, from bottom to top, type II migration timescales 
for a ss = 1, 0.1, 0.01, 10 -3 , 10 -4 , 10 -5 respectively. 
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Figure 11. The minimum initial eccentricit y required for ec- 
centr icity evolution via planet-disk interaction iGoldreich fc Saril 
2003). Heavy line is for a Jupiter mass planet and light line is for 
a 10 Jupiter mass planet. Obviously, eccentricity growth is more 
likely for a planet in the dead zone. 
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Table 1. Recombination coefficients in rate equations. 



symbols 


recombination rate coefficients [cm 3 s J ] 
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APPENDIX A: EFFECT OF GRAINS ON DISK 
IONIZATION BALANCE 

To calculate an electron fraction by including various recom- 
bination sources, we write the rate equations for electrons, 
molecular ions, neutral and singly charged grains as follows: 



— — — £n 

dt s 



dn„ 



dt 



dn G 



P n m +n e — p r n M +n e — Pi n G ne 
-Pi n G +n e (Al) 
= £w„ — p n e n m + — Pt nMn m + — Pi nGn m + 



-Ps n c 



(A2) 



—rr = P2 n a +n e + p 5 n G -n m + + p 7 n G -n M + 



dt 
dn n 



+2 Pi n G +n G - — Pi n c n G — Pi n m+ n G 

—Pa n M+ n G (A3) 

= Pi n G n m + + pe n,Gn M + — Pi n e n G + 

-P 3 n G -n G+ (A4) 



dt 



= Pi n G n e - p 3 n G+ n G - - p 5 n m+ n c 



—Pi riM+n G - 
n e — n m+ + riM+ + n G+ 



(A5) 
(A6) 



Here, n e , n m , um, no, and n„ are the density of electrons, 
molecules, metals, grains, and neutral hydrogens, while the 
superscript of + or — indicate the elements are positively 
or negatively charged respectively. Also, £ is the total ion- 
ization rate /3s are recombination rate coefficients (see Ta- 
ble for details). We solve these differential equations to- 
gether with charge conservation by using a semi-implicit ex- 
trapolation method (a genera lized Bulirsch-Stoer method) 
(stifbs.f90 in lPress et alll 996l . This calculation takes about 
a few hours to run on the average workstation. 

Table^shows all recombination coefficients appeared in 
above equations, where fcs is the Boltzmann constant, T is 
the disk temperature, a is a grain radius, and q is the charge. 
Also, m e , m m , mja, and m G are the mass of electrons, 
molecules, metals, and grains respectively. The recombina- 
tion ra te coefficients P, p r , p t are taken from lFromang et alJ 
( 2002) , and the rest are obtained by following the method 
bv lSano et aflfeOOCh . 



e > 



(B2) 



where w is a gap width and M* is the stellar mass (SG04). 
Combining the second and third points, we can write 



M d r p M d 



(B3) 



where Md is a disk mass. The left inequality arises because 
a gap has to open slower than the eccentricity evolution 
timescale, while the right inequality comes in so that the 
significant decay won't happen. 

Writing these as the minimum and maximum ratios of 
a gap width and a disk radius, we can rewrite first two con- 
ditions as follows: 



e > ei 



2/3 



> 



Mp/M, 
e > ei 

M v 



5/3 / X 4/3 
W \ I tip 



e2,j 



^1 ,min 



(B4) 



(B5) 



Thus, the eccentricity must be at least e > 
max [ei iin i n , e2 jm i n ] to be excited, which is Equation Q. 



APPENDIX B: MINIMUM INITIAL 
ECCENTRICITIES 

In §5, we argued that there is the minimum initial eccentric- 
ity required to e nhance a planetary ecce ntricity via planet- 
disk interaction. ISari fc Goldreicbl i2004f) suggested that for 
the eccentricity to grow, (1) the corotation resonances must 
saturate, and (2) eccentricity must grow faster than a gap 
opens. Also, they claimed that eccentricity does not signifi- 
cantly decay as long as (3) corotation resonances are more 
than 5 % saturated. 

For the first point, two equations must be satisfied: 
tsat < t V i S ; which may be written as 

2/3 / \ 5/3 / v 4/3 

«>l^fe) ft) (B1) 



as well as t sa t < tgap', which becomes 



